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Abstract 

We consider a model describing the behavior of a mixture of two incompressible fluids 
with the same density in isothermal conditions. The model consists of three balance 
equations: continuity equation, Navier-Stokes equation for the mean velocity of the 
mixture, and diffusion equation (Cahn-Hilliard equation). We assume that the chemical 
potential depends upon the velocity of the mixture in such a way that an increase of 
the velocity improves the miscibility of the mixture. We examine the thermodynamic 
consistence of the model which leads to the introduction of an additional constitutive 
force in the motion equation. Then, we prove existence and uniqueness of the solution 
of the resulting differential problem. 

AMS Classification: 35Q35, 76T05. 
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1 Introduction 

We consider a model describing the isothermal motion of a mixture of two incompressible 
fluids following the diffusional approach to binary mixtures. This goes back to model H 
in the classification by Hohenberg and Halperin [12], which consists of a Cahn-Hilliard 
diffusion model coupled with a fluid motion. This kind of approach is extensively discussed 
by Lowengrub and Truskinovsky |14j . Basically, a binary mixture is described in terms of 
a macroscopic velocity field representing the mean velocity at each spatial point, the total- 
density field of the mixture, and an order field which describes the actual composition of the 
mixture at each point. The model consists of three balance equations: continuity equation, 
(total) momentum balance (Navier-Stokes equation), and diffusion equation (Cahn-Hilliard 
equation). The diffusional model can be considered as an approximation of the classical 
theory of mixture (based on two continuity equations and two momentum balance equations, 
one for each component) when the momenta and the kinetic energies of the relative motion 
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can be neglected. 

The coupling of the motion equation with the Cahn-Hilliard equation has a trivial part due 
to the presence of the material derivative (rather than the partial time derivative) of the 
order parameter in the Cahn-Hilliard equation and a less trivial coupling arising from the 
dependence of stress tensor upon the gradient of the order parameter. The presence of such 
a coupling has been derived by Gurtin et al. on the basis of classical continuum mechanics 
arguments 

In this paper we discuss a variant of this model in which a dependence of the chemical 
potential upon the velocity of the mixture is introduced. We add the velocity-dependent 
term in the local part of the chemical potential, that is the part independent from the 
gradients of the order parameter. The effect of the velocity can be assimilated to an increase 
of the temperature (which is a fixed parameter in the isothermal model we are considering), 
that is it reduces the miscibility gap. 

In section [5] we review the classical analysis of the thermodynamic consistence which 
displays the need of an additional constitutive force term in the motion equation. 
The subsequent sections are devoted to the proof of existence and uniqueness of the solution 
of the resulting differential problem. Our mathematical study concentrates on the fully 
incompressible situation, that is the case of a binary mixture of two incompressible fluids 
which also have the same density. Clearly, this is an exceptional case from an empirical point 
of view, but could be an acceptable approximation for a broader class of real situations. 

The coupling of the Navier-Stokes equations with the Cahn-Hilliard equation has been ex- 
tensively studied in the literature. Among the first results concerning existence, uniqueness 
and asymptotic behaviour of the solutions, we recall the papers [3J and [16]. More recently, 
Abels [1] proves well-posedness and examines long-time behaviour of the Cahn-Hilliard- 
Navier-Stokes system involving a class of singular free energies, (including the logarithmic 
free energy) which guarantees the boundedness of the order parameter. Further results 
about the asymptotic behaviour of the solutions and the existence of global and exponential 
attractors for the coupled system are shown in [9]. 

The main result of our paper is the proof of existence and uniqueness of the solutions 
of Cahn-Hilliard-Navier-Stokes equations where a new non-linear term is present due to 
the velocity-dependence of the chemical potential. Moreover, we add a viscous term in the 
definition of the chemical potential (|10j) which turns out to be crucial for our purpose. 
The functional formulation of our problem is given in section [31 In section [H with the 
same technique used in [8], we introduce a family of approximating problems, by adding 
to the Cahn-Hilliard equation a perturbative term proportional to the time derivative of 
the chemical potential. By means of a fixed point argument, we establish the existence of 
solutions of the approximated problems. Finally, in section [5] we prove well-posedness of the 
original problem letting the perturbative term tend to zero. 

2 Model equations and t her mo dynamical consistence 

Let us consider a mixture of two partially miscible fluids; we will use a binary index i = 1, 2 
to make reference to each of them. Every spatial volume element dV will in general contain 
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a mass portion drrii of the i-th fluid; we indicate with pi the apparent density of each fluid: 

drrii 

The adjective 'apparent' is used to emphasize that we are considering the ratio of each mass 
fraction over the total volume element dV, rather than over its own fractional volume dVi. 
Each component, which can be compressible or incompressible, is characterized by its own 
density p® = drm/dVi in standard conditions of temperature and pressure. Of course, the 
total density is the sum p = pi + pi. We also define an order 'parameter measuring the 
degree of phase separation as 

dm,! - dm 2 pi - p 2 , , , , 

f = l r~5 — = e -!> 1 • 

dm i + dnri2 p 

As we are considering mutually non-transforming chemical species, the first general balance 
laws we have to impose are the mass conservations of each component; so, if Vj is the velocity 
of the i-component, we demand 

^ + V-( ft vO = > i = 1,2. (1) 
By defining the mean velocity of the mixture as 

V = p~ X {pl\l + /9 2 V 2 ), 

(so that pv is the total momentum density), a global mass continuity law follows 

^+V-(pv) = 0. (2) 

The model we are going to study describes an incompressibile mixture of two fluids, so the 
total density p = p\ + p 2 is constant and the continuity equation ^ reduces to 

V-v = 0. 

This means that each fluid component is an incompressible fluid and also that each compo- 
nent has the same constant density p: 

drrii = pdVi. 

From equations ([T]) it also follows, after a little calculation, the following equation for ip: 



(v 2 - Vi) 



We will use the usual notation for the material time derivative related to the mean velocity 
field v (barycentric material derivative): 
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for every spatial field f(x,t). Therefore the equation for ip assumes the form 

pip = V • J. (3) 

This is the usual equation for a conserved quantity with respect the gross motion of the fluid 
defined by the mean velocity v. 

We introduce now a basic physical hypothesis which characterizes the diffusional approach 
to binary fluids f [111 114] ). Accordingly, we will describe the dynamics by using a balance 
law for the total momentum density pv of the mixture, while the effects of the relative 
motion will be described only through the balance law ([3]) for the scalar order field ip, not 
retaining the motion equation for the relative momentum P2V2 — pl^i- In other words, the 
fundamental fields of the model will be p,v,ip rather than px,P2i Vi,V2, and the current J 
will be given a constitutive law in terms of p,v,(p (and their gradients). Physically, this 
amounts to neglect the kinetic energies and the momenta of the constituents relative to the 
mean motion, only retaining the information of the relative mass flux. The resulting model 
can be considered as the model of a single fluid with an internal defined by the conserved 
field ip. 

So in this paper we are going to consider a model characterized by the three balance 
equations: 

V-v = 0, 

pv = V • T + pd + pi , (4) 
pip — V • J. 

Here T is the stress tensor, d a constitutive body force and f a possible external body force. 
The (unusual) constitutive force d is required if we want to construct a thcrmodynamically 
consistent model in which a dependence of J on v is admitted, as we are going to show. 

We begin considering the diffusion equation. We will consider a Cahn-Hilliard similar 
model, in which the current is expressed as 

J = 7 V M , (5) 

where p is the chemical potential. In the classical Cahn-Hilliard model ([2J [5]) p is a non 
local function of the order field ip which takes the form 

P = MiocM - V ' 

where h = Vip. Here we want to consider a generalization which possibly includes a de- 
pendence on velocity v. Also, for future utility in the mathematical study, we will add a 
dissipative contribution proportional to ip. So we assume 

P = A*ioc(v, ip) - «i(v?)V • [K 2 (ip)h] + /3ip. 

We associate a suitable balance of powers to the diffusion equation, which is needed to 
write the first law of thermodynamics later on. This is obtained by multiplying the diffusion 
equation by the chemical potential p 

pipp = /iV ■ J. (6) 
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The central issue is to recognize in this equality an internal and an external power. The 
choice will be influential in satisfying the second law of thermodynamics. In particular, 
as there are non local (gradient) contributions in fx, it would create difficulties to refer 
completely the term pipp to the internal power. 
We will rewrite equation in the form 

ptfUioc + j3p<p 2 + ■ V(p«i<p) + J • V/i = V • (pJ + pniK 2 iph) 

and we ascribe the left hand side to the internal power "P^H Letting L = Vv (in components: 
Lij = dvi/dxj), we observe that 

V<p = h + L T h. 

We use this identity and ([S]) to write the internal power in the more suited form 
K = [pMloc + « 2 h ■ V{pKi)\<p + Pptp 2 + pn lK2 h • (h + L T h) + 7 |Vp| 2 . 

Next we consider the momentum balance equation. The stress tensor and the constitutive 
force are given by 

T = T(D,p,h), d = d(v, 

where D := 4(L + L T ) is the symmetric part of the velocity gradient. 

The balance of powers is obtained by multiplying both members of (j4])2 by v: 

ip(v 2 )- = [V ■ (Tv) + pf • v] — [T : Vv — pd ■ v]. 

The term d • v from the constitutive body force will contribute to the internal power. In 
particular, the internal mechanical power V l m is defined by 

Vl = T : D - pd • v. 

This identification is based on the assumption that stress tensor is a function depending 
only upon the first gradients of the field, namely T = T(D, ip, h); otherwise, as it happens 
for /i, which is dependent upon Vh, it would be appropriate to refer part of the contribution 
T : D to the external power. 

As we are considering an isothermal model, we use the dissipation inequality 

(4-vi- vl < o, 

where ip the free energy, as the proper version of the second law of thermodynamics. 

If ij) = ip(x), where x is the list of the variables v, ip, ip and all their gradients appearing in 

the constitutive equations, the dissipative inequality is written as 

P ^2 ^XiXi + \pipif - PP\ oc - «2h • V(pKi)] (p - f3p(f> 2 + [p^h - p«iK 2 h] • h 
-[T + pKiK 2 h(g)h] : D + pd- v - 7 |Vp| 2 < 0, 

1 There is an alternate approach to the issue of the energy balance associated to the diffusion equation due 
to Gurtin. According to this author, the diffusion equation is associated to the internal power pipfi + J ■ V^t, 
which could be considered as a definition of the chemical potential, which is treated as an independent field. 
Moreover the diffusion equation is accompained with an independent balance equation (microforce balance) 
which also brings a contribution to the total amount of power. 
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having used h • (L T h) = D : (h ® h). It easily follows that ip does not depend upon any of 
the variables Xi ¥^ ¥■> h, that is 

V> = ^,h). (7) 

For any given tp, h, it ispossible to find processes with tp = 0, v = 0, D = 0, V/x = but 
otherwise with h arbitrary^ This implies 

iph = KiK 2 h. (8) 

By choosing appropriately h we can make V/x = with tp, tp, h, v, D arbitrary, so the 
following inequality holds 

[pi/V _ /°Mloc - K 2 h. ■ V(pKi)] tp ~ fiptp 2 - [T(D, ip, h) + pKi^h. ® h] : D 
+pd(v, tp, tp) • v < 0. (9) 

In the same way, we can make V p ^ and tp = 0, v = 0, D = 0, so 7 > 0. 

Inequality ([9]) implies that j3 > (considering processes with v = 0, D = and tp arbitrary); 

letting tp — and v = we have to impose 

[f (D, tp, h) + pKi« 2 h <g) h] : D > 0. 

Because of the incompressibility, on one hand the pressure (that is the trace part of T) 
is not a constitutively determined quantity (it is kinematically determined), on the other 
hand the trace part of D identically vanishes 

Tr(D) = V-v = 0. 

So, putting for brevity, 

Q = T + pKiK 2 h(g)h, 

the inequality Q : D > is equivalent to Q : D > where the tilde on a tensor indicates 
its deviatoric part: D = D — -|Tr(D)l and similarly for Q. So we assume Q = 2vT> with 
v = 0(T>, tp, h) > 0, that is 

T + pKi/c 2 h (g) h = 2z/D - pi. (10) 

where p — — ^Tr(Q) is the indetermined component of the pressure. 
We are now left with the reduced inequality (for D = 0) 

\fnp<p(cp, h) - pp,i oc (tp, v) - K 2 h ■ V(pKi)] tp - fiptp 2 + pd(v, tp, tp) ■ v < 0. (11) 

From that we see that, if there is a non trivial dependence of p,\ oc on v, the presence of 
the constitutive force d is necessary. The following particular choices are made to satisfy 
(HU): 

d = 8{w,tp)tp, (12) 
Ki = constant, ip v h — 0, (13) 
Woc(p, <p, v) - 6 ■ v = tp v (tp) (14) 



2 In fact V^t = V(/^i oc + f}tp) + V[fC2V ■ («ih)], so one can make V^t = by suitably choosing VVh for 
the given set of conditions. 
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Conditions 0, ©, ([TD ]) .([T2 j) -([H )l with /3 > 0,7 > 0, ensure that the dissipation inequality 
is satisfied. 

The specific feature of the model we are going to study (compared with those by Gurtin 
and Truskinovsky) is a velocity dependence of the chemical potential. We adopt the usual 
form of the free energy function ip(ip, h) used in the Cahn-Hilliard model of diffusion: 

^ = ^^ + u-G^) + H{ v ) 1 G(<p) EE itp 2 , i^EE-^ 4 , 

where u G [— l,+oo[ is a temperature-dependent parameter, typically 




The ^-dependent part of ip is such that for u > it has a unique minimum at (p = 0, 
while for —1 < u < has two minima ±(p with (p s]0, 1]. It is known ([5]) that the unique 
minimum in the potential corresponds to the situation without miscibility gap, while in the 
regime with two minima there is a miscibility gap. 

In this paper we assume the following form for the local part of the chemical potential 

AiocO, v) = i> 9 {ip) + Xv 2 G'(ip). 

We remark that the effect of velocity can be assimilated with an increase of temperature, 
that is a restriction of the miscibility gap. So <5v = Av 2 G'(ip) and we obtain the constitutive 
force 

d = Av G'{y)<p = \vG(cp). 

We sum up the system of equation in that case, putting everywhere p — 1 and taking 
7, v, k — K1K2 > constant: 

V-v = 0, (15) 

v = -Vp + ^Av - kV • (V<p ® Vip) + Xiptpv + L (16) 

ip = 7 A M , (17) 

where 

/i = -nA(p + ip 3 + (u + \v 2 )tp + fi<p. (18) 

3 Notation and functional settings 

In order to obtain a precise formulation of the problem, we introduce here some notation 
and recall the main inequalities used in the sequel. We assume that the domain il occupied 
by the system is a bounded subset of R 2 , with smooth boundary dil. 

For each p > 1 and s <E K, we denote by L P (Q) and H s (fl) the Lebesgue and Sobolev 
spaces of real valued or vector valued functions, according to the context. Let || • \\ p and 
|| • \\hs be the standard norms of L p (f2) and fP(f2), respectively. In particular || • || stands 
for the L 2 (f2)-norm. The space Hq(£1) is the closure of C°° functions with compact support 
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with respect to the norm || • Wh 1 - Moreover, is the dual space of endowed 

with the standard norm 

\\w\\ {H iy =sup{\{w,u)\ :u€H 1 (Q.),\\u\\ H i < 1}, (19) 

where (•, •) denotes the pairing. 
We define 

H\n) = {weH 1 ^) : Vu;-n| a n = 0}, 

H 2 (n) = {weH 2 (n) : Vwn\ dQ = 0}. 

For vector valued functions, we introduce the functional spaces used in the framework of 
Navier-Stokes equations (|18j) 

Hi v (n) = {weif^iV-w^}, 

L 2 W {Q) = {weL 2 (r!):V-w = 0,vn| 9n = 0}. 

Finally, for any T > we define 

X v = L 2 {0,T;H 3 {Q)nH 2 {Q))r\H 1 {0,T;H 1 (n))r\H 2 (0,T;H 1 {ny) 

X v = L 2 (0,T; J ff 2 (r»)ni/ d 1 iv (fi))nH 1 (0,T;L2 iv («)) 

X^ = L 2 (0,T-H 2 (Q)) 

Xt — Xip x X v x X^ 

endowed with their usual norms | • \\x v , \\ ■ \\x v , \\ ■ || x„ and 

ll(^v,/i)||^ T = |M|^ + ||v||i v + |H|^. 

Here and henceforth we denote by C any positive constant depending only on the domain 
n which is allowed to vary even in the same formula. Further dependencies will be specified. 
Since del 2 the Sobolev embedding theorem implies ([2]) 

\\w\\ p < C\\w\\ H i, 2<p<oo, w£H\Q) (20) 

and the following interpolation inequalities hold as a consequence of the Gagliardo-Nirenberg 
inequality [71 [T5] : 

HI* < c-HHIIkll^, (21) 
Hie < c|M| 4 /3|M|^ 8 . (22) 

If w € Hq(Q) or w € iJ 1 (f2) and J Q wdx — 0, the Poincare inequality provides ([B]) 

\\w\\ < C\\Vw\\ . 
From the Agmon inequality ( jTTJ p.52]), we deduce that 

\\w\\ x <C\\w\\ H 2, w£H 2 (n). (23) 
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Furthermore, for every v G H (O), u, w G H 2 (fl) the following interpolation inequalities 

< qiwIlffiHltfllfla (24) 

||uw||^2 < C||m|| h 2 \\ w \\h 2 (25) 

hold as a straightforward consequence of (|2Tfl) and (|23p . 

In addition, if 10 G H 2 (fl) n fl^(n) or it; G ff 2 (ft), then ([H Thm. 5.1, pag. 149]) 

\\w\\„2 <C(\\w\\ + \\Aw\\). 

As a consequence, for every w G H 3 (£l) n H 2 (£l), we have 

\\w\\ H 3 < C(\\w\\ H i + \\Aw\\ H i) < C(\\w\\ + \\VAw\\). (26) 

For vector- valued functions we define the orthogonal projector P : L 2 (Q) — > L^ iv (0) and 
the operator A defined as 

Aw = -PAw, w G H 2 {Q) n i?div(°)- 

It is worth noting that for any w G H 2 (£l) n H^ iv (Sl) the following estimate holds: 

j|wj| H 2 < (7(||w|| + ||Aw||). 

For later use, we will also need the following result, whose proof is given in [51 Thm 4, 
pag. 288]. 

Theorem 3.1 Suppose w G L 2 (0, T; H m+2 {Vt)) and w t G L 2 (0, T; H m (il)) where m is a 
nonnegative integer. Then w G C([0,T]; H m+1 (fl)) and 

mS^\\w(t)\\ Hm + i < C(|MU 2 (0,T;,ff'"+2(fi)) + \\wt\\L*(0,T;H<"((i))), 

the constant C depending only on T, f2, and to. 

Finally, for reader's convenience, we recall Young's inequality. Let 1 < p,q < oo, such 
that 1 + 1 = 1. Then, 

ah <r]a p + C{r))b q , (a, b > 0, r) > 0). (27) 

The functional formulation of problem (fT5 l) - ([T8| is the following: 
Problem (P). To find a triplet ((p, v, fi) G Xt such that 

/8ip t - nAip - fi + (3v ■ \7tp + <p 3 + utp + Xipv 2 = (28) 
v t + zMv + P [kV • (V</3 ® Vy>) - A<p(pi + v • V<^)v + (Vv)v] = Pf (29) 
-7A^i + ipt + v ■ = (30) 

a.e. in Q x (0, T), and 

(p(x,0) = (po(x), v(x,0) = Vo(x), a.e. x G fi. 
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Notice that equation ([25]) which provides the definition of the chemical potential, is inter- 
preted as a parabolic equation governing the evolution of ip, whereas equation (|30[) is an 
elliptic equation for the unknown fi. Accordingly, the condition /3 > will be crucial in the 
following analysis to prove well-posedness of the system. 

Existence of solutions to problem (P) is proved by introducing a suitable family (P £ ) of 
approximating problems, where e is a small parameter such that < e < 1. In Section U we 
prove existence of solutions to (P e ). Then, by letting e — > 0, we obtain the existence result 
for the solutions of problem (P). 

4 Approximating problem 

We construct a family of approximating problems of (P), by adding the term £/x t to equation 
([30]) . Accordingly, we introduce the functional space 

Xj- = X<p x X v x A E , 

where 

X £ = L 2 (0,T;H 2 (n)) HP^O,! 1 ;! 2 ^)). 

Problem (P E ). To find a triplet (y> E ,v E ,/x e ) € Xf, such that 

/3y; E - kA(^ £ - fi e + j3w £ ■ \7ip £ + {tp £ f + upf + \ip £ (v £ ) 2 = (31) 
v t e + vAv £ + P[kV • (Vip £ ® Vy> e ) - A</(<^ + v £ ■ V</)v e + (Vv £ )v £ ] = Pf (32) 
e/4 - jAfi £ + (p £ + v £ - Vip £ = (33) 

a.e. in fl x (0,T), and 

<p(x, 0) = ip (x), v(x, 0) = v (a;), [i(x, 0) = Mo(z)> a -e- £ G 0, (34) 

where /^o is the solution to equation 

-/3jAno + /i = -kAi^o + ucpo + A^oVg. 

From the standard theory of linear parabolic equations (see e.g. [HI [T3J [18] ) , we deduce 
the following auxiliary result. 

Lemma 4.1 (Existence of solutions of problem (Pl)) Let $ G L 2 (0, T; P 1 (fi)) 
nffHO/;]? 1 ^'), T G L 2 (0,T;4(I!)), A G L 2 (0, T; L 2 (fl)), ^ G P 2 (0) ; v G ffi v (fi), 
/Lto G P x (ri). Then, there exists a unique solution ((/?, v, /x) G Af. of the linear problem (Pl) 

P<p t - = $ (35) 

v t + uAv = T (36) 
efit-jAfj, = A (37) 

witt ifte initial conditions In particular, tp G C([0, T]; H 2 {tt)), v G C([0, T]; P] iv (0)) ; 
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Lemma 4.2 Let C0,w,C) G Xf,, f G £ (0, T; L 2 (fl)). Then, the functions $, T, A de/med 



fl.s 



$ = C - /3w • V?A - V- 3 - - A-0w 2 (38) 
T = 7>[f - kV • (VV'O Vr/0 + A^C^i + w • V^)w- (Vw)w] (39) 
A = -i/j t - w • VV> (40) 

safe/y $ G L 2 (0,T;i/ 1 (^))ni/ 1 (0,r; J ff 1 (^)'), T G £ 2 (0, T; Z^ iv (0)), A G L 2 (0, T; L 2 (f!)). 
Proof. From definition (j3"8")l and inequalities (IM1) and (j2"5)l , it follows that 

11*11//! < C(||C||Hi+||w-V^||Hi + ||^ 3 ||Hi + ll^llHi + II^W 2 || ffl ) 

< C(\\C\\m + \W\ H ,W\\ H , + IIVII^HVII^ + UWrn 

+||V'IIh2|I w IIh i II w IIh 2 )- (41) 

The assumption (ip,w, Q G ATf. and Theorem 13.11 guarantee that ip G C(0,T; H 2 (il)), 
w G C(0,T;^ iv (Q)), C G C(0, T; H Accordingly, Young's inequality yields 

T 

\m 2 H idt < co. 

In order to prove that $ G i? 1 (0, T; H 1 ^)'), we differentiate equation with respect 
to t and we evaluate ||<I> t ||(#i)/. In view of ([TT))) and the regularity of the functions tp, w, £, 
we have 

ll^lkffi)' < C(||Ct|| + ||w t ||||vv||fli + ||w|| ffl ||w^|| + UWhWM + \\M 

k\\W\m + \\1>\\& l|w||Hi||w t ||). (42) 



In addition, Young's inequality leads to 

WQtWfaydt < c£[\\Ct\\ 2 + (1 + Hwllir, + W\UWt\\m 

+(11^11^1^11^ + ||V||^)l|w t || 2 ]di< 00. 

Now we consider equation (p?9")) . We get 

||T|| < C(||f|| + U\\ H 2U\\ H 3 + \\iP\\ m U t \\m\H\m + \\1>\\%* + M\hAW\h*), 

which implies T G L 2 (0, T; L^ iv (fi)). 
Similarly, from (l4"U|) we deduce 

||A|| < HVtll + ||w|| H i||^|| K 2, 

so that A G L 2 (0,T:L 2 (n)). □ 

Theorem 4.1 Suppose that tp G H 2 (fi), v G #div( ft )> Mo G H 1 ^), f G L 2 (0 , T ; L 2 {Q)) . 
Then, problem (P e ) admits a unique local solution for a sufficiently small time t G (0,T). 
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Proof. For any r > 0, we define 

S:X F - -> X% 

i->- (^,v,/i), 

where (y>, v, /x) is the unique solution of problem (Pl) and $, T, A are defined by (|38l) - (|40]) . 
Thanks to Lemmas 14. 1114. 21 S is well defined. Furthermore, we consider the bounded subset 
B T C X^ that consists of the functions (ijj, w, £) satifying the following conditions: 

h + \H\ 2 H 2 + WCWUdt < ni , (43) 

Mm+\\™ t \\ 2 + \\Ct\\ 2 }dt<n 2 , (44) 

Htt\\fmy dt < n 3, (45) 

where «i,n2,ri3 are positive constants which will be specified in the sequel. In particular, 
in view of Theorem 13. 1[ there exist suitable constants C(rii,n 2 ), C(n 2 ,n 3 ) such that 

rnaxOI^H^ + \\w\\ 2 m + HCH^x] < C(ni,n 2 ), (46) 

[0,r] 

maxll^H 2 <C(n 2 ,n 3 ). (47) 

[0,r] 

The proof is divided into two steps, 
i. 5 maps B T in itself. 

Throughout this proof we denote by C a generic positive constant which is allowed to depend 
also on e. Let us consider equation (f35j) . By multiplying in L 2 (f2) by ip and integrating by 
parts, we obtain 

H||^ + K ||v^|| 2 <^||$|| 2 + C||^|| 2 , 

where r\ is a (small) positive constant which will be chosen later. From definition (|38p of <3>, 
the Sobolev embedding theorem (|20p and relations (|46l) - (jT7)) . we deduce the inequality 

||*|| < cflKII + ||w|| ffl ||vv|| ffl + HVIIhi + IMI + MIjpI|w|&i) < c(m,n 2 ). 

Hence, an application of Gronwall's inequality leads to the estimate 



M| 2 < ^ 



\vo\\ 2 + 2 fl in 2 * 



<e CT [\\<p f + C( ni ,n 2 )r,T\ 



Choosing m > 2||< / 9o|| 2 and r\ and r small enough, we infer that ||</?|| < n\. 

Now we multiply (|35|) in L 2 (Q) by (</?t — Aipt) and we integrate by parts. By taking 
Young's inequality into account, we obtain 

~l\\^\\ 2 + HV^II 2 ] + ^\\ Vt f H1 < C||$|| 2 K1 . (48) 
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From (pfTj) . we deduce that 

11$ 11^ < C(\\cf H > + llw • v^|| 2 ffl + HVII^IIV'II^ + IHI 2 ffl + ll^w 2 ||^)- 

In particular, in view of (|20l) . (HH, (P251) . and Young's inequality, we have 
llw-V^Uffi < C(||w V^|| 2 + ||(Vw)W|| 2 + ||(VVV)w|| 2 ) 

< c(||w|| 2 ffl ||^||^ + livwii 2 !!^!^ + livwiffllwii^) 

< CdlwH^I^H^ + HV-ll^llwll^ + ||Vl| 2 ^l|w||^) 

+?7l|w||^ 2 +??||V'IIh3, 

for any rj > 0. By repeating similar arguments, we infer that 

H™ 2 \\m < C(||^w 2 || 2 + ||w 2 W|| 2 + ||^(Vw)w|| 2 ) 

< Cd^ll^llwH^ + HVVll^llwll^ + IIVII^IIVwIIIIIwII 2 ^) 

< cmUM^m + ||^|| 4 ffl ||w|| 6 H1 ) + v\Mh- 

Collecting (|48|) - (l5TT) and taking (I43)) , p6)) . into account, we prove the estimate 

|[IIa^ii 2 + iiv^ii 2 ] + ^£y t \\ 2 m dt 

< ^[HA^oll 2 + HV^oll 2 ] + 3r?ni + C{n u n 2 )r. 

Thus, taking n 2 greater enough according to HyoH^ and r, r\ small enough, we have 

/ \\tpt\\ 2 H idt <n 2 . 
Jo 



From (J35J) it follows that 

it\\A<p\\ m <P\\vt\\w 
Hence, on account of and we deduce that 

IMI#3 < Til, 

Jo 

with r small enough and rii > n 2 . 

Finally, we observe that by differentiating (|35p with respect to time, we obtain 

(itpu - nAtpt = $ f , 

which implies 



\Wtt\\i my dt<C I [11^11^+11^11(^1) 
Jo 

Owing to (|32 ]) -(|35 ]) . we get 

WvttWfH 1 )'^ ^ C(n!,n 2 ) + Cr. 



]dt. 
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Choosing r small enough and n 3 > C(rii,7i2), we deduce that 

\\<Ptt\\\ H iydt < n 3 . 

Now we multiply (|36[) by (v + Av) and we integrate over f2, thus obtaining 
i|lHI^ + ^l|Vvf + |||Avf< C ||Tf. 
The definition (|3l?l) of T implies the following inequality 

l|T|| 2 < qi|f|| 2 + ||vv^||2||v^||l + ||^|| 2 ffl ||v. t |lll|w||^ 
+II^IIh 1 ||w||^||VII^ + I|Vw||2||w|| 2 ] 

< vWWh* + vUtf Ht +^||w|| 2 j2 +C||f || 2 +C(n 1 ,n 2 ,n 3 ), 
for any 77 > 0. An integration of (I5^t over (0, r) leads to 

Jo 

< II v o|Ihi + »7(2ni + n 2 ) + C(ni,n 2 , n 3 )r + C\\f || 2 2 (0iT;i 2). 
We take m, n.2 large enough and 77, r small enough. Accordingly, 

HvH^dt < C [ T (\W\\ 2 + \\Avf)dt<n u 
Jo 



and a comparison with (|36[) yields 

|v t || 2 ^ < n 2 . 







We multiply (|37| by (zx — A/x) and we integrate over fi, thus obtaining 

f |ll A + JIAmII 2 + 7IIVHI 2 < C(||/x|| 2 + ||A|| 2 ). 
In view of (gDJ), (gBJ) and (gTJ), we infer that 

II A|| 2 < ||V^|| 2 + C||w|| 2 ffl ||^|| 2 ff2 <+C( ni ,n 2 ,n 3 ). 
As a consequence, by applying Gronwall's inequality to (|54[) we obtain 

§IHI!p + \f \\^\\ 2 dt < Ce CT [HmoIIIi + C(n 1>na , «3)r] . 

By choosing n\ large enough and r small enough, we deduce that 

f \\n\\ 2 H 2dt<C ( (II/jIP + IIA/xII 2 )^^^!. 
Jo JO 
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A comparison with (|57|) yields 

£ 2 || Mt || 2 <2( 7 2 ||A / i|| 2 + ||A|| 2 )< 2 7 2 ||A M || 2 + CK,n 2 ,n 3 ). 
Accordingly, we have 

||^f|| 2 ^ < n 2 , 

provided that ni is large enough and r is small enough. 

2. S is a contraction in B T if r is small enough. Let {ipxi Wi, Ci), (^2; w 2, C2) £ -St and 
(<^i, vi, /ii), ((^2, V2, H2) be the corresponding two solutions of the linear problem (Pl) with 
the same initial data. We denote by ip — ipi — 1P2, w = wi — W2, C = Ci — C2 an d prove that 
S : B T — > B T is a contraction mapping with respect to the metric induced by the norm 

HI (V, W, C)IH 2 = IMl£a(0,T;ff») + \\Mh(p,T-,&) + \\Mo{[0,t];L') + II W H h(0,r;H 2 ) 
+ II w *IIl 2 (0, t;L 2 ) + IICllc([0,i-];.ff 1 )- 

It is worth noting that X T is a closed metric space with respect to the metric 1 1| • 1 1| . Therefore, 
our aim consists in proving that 

|||5(Vi,w 1 ,Ci)-5(^ 2 ,w 2 ,C 2 )||| <i T |||(^w,C)|||, (55) 

with < L T < 1. From (|35|) and (f38j) . we deduce that ip = ipi — ip 2 satisfies the following 
equation: 

/3cp t - nAtp = $1 - $ 2 , (56) 

where 

$1 - $2 = C - /?( w i • W> + w • Vip 2 ) - ip{^l + ipi^2 + - uip 
— X[tjjw1 — ^(wi + w 2 )] • w. 

Let us multiply (|56[) by <p and integrate over fi. By means of (|20[) and Young's inequality, 
we obtain 

^||MI 2 + K ||v^|| 2 < c|M| 2 + r, [IKII 2 + (ll^ill 4 ^ + ll^ilMV* + UAm 

+l|wi||iri + l)U\\m + (11-02 1|^ 2 + II^IIhiHw! +W2|| 2 ^ 1 )||w|| 2 /1 + ||w 1 || 2 fl ||V>|| 2 ^] 
for every 77 > 0. The assumption (ipi, Wi, Ci), (^2, W2, C2) G £? T guarantees that 

^IMI 2 + K ||v^|| 2 < qi^ii 2 + c v [HOI 2 + M 2 H , + \H\U . 

Hence Gronwall's inequality leads to the estimate 

IMI 2 < C V e CT [ T [HCII 2 + + \\wf m ] dt < CVe CT |H(^w,C)||| 2 . 

Jo 
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We multiply (|55|) in L 2 (f2) by (ip t — A(pt)- An integration by parts and Young's inequality 
yield 



iivvn+|iiv,n 2 s , < cii#i 

From the definition of $1 — $2 and inequalities (|24|) . (|25|) . we have 

|$i - $ 2 ||ffi < c [HCIIffi + ||wi • vvil^i + ||w • v^ll^i 



$a||ln. 



$2 

4 



J 2||ffi ^ ^ Lll^llH 

-(!l^ill 4 ff2 + l^ill^HV^II^ + ||^|| 4 ff2 + 



+IIV' w illfl-i + IIV>2|Ih=II( w i + w a) • w llffi] • 

Owing to Sobolev embedding theorem and accounting for (f^Uj) . (p?Tj) and (HrJl) . we obtain 

|| (wi +w 2 ) • w||^i < C(||wi +w 2 ||^i||w||^ 1 
+ j|V(wi + w 2 )||i||w|| 2 ?1 + ||wi+w 2 ||^i||Vw||J) 
<C(\\w\\ 2 m 



2 

< C(l 



||wi + w 2 || H 2||w|| H1 + ||w|| H i||w|| H 2) 

|Wl|| H 2 + ||w 2 ||jj-2)|| w|||^i + 7]l||w||^2, 



with rjx > 0. Moreover, proceeding as in the second inequalities of ([5Hf and (|5Tj) we can 
prove the estimates 



lWH 1 



|wi • Vlp\\ H l 



< 



< 



h-2 + HwiHja-2 HVllja-i), 



c{U\\l 
c(i + 1^11^)11^1 



2 

H 2 



2 

H3, 



llw-V^H^ < C(l + U2\\m)\\M\ 2 H ^ +ml|w||ff2, 
where f}x,r}2 are suitable positive constants. Collecting the previous results, we get 

f |[(I|A^I| 2 + ||V^|| 2 ) + ^||^|| 2 ffl < 2r 7l ||w||| f2 + m \\n% 3 

+ C [llC|lc([0,T];Hi) + C 1 + ll W l||ff 2 )IIV'llc([0,T];^) 

+ (1 + l|wi|| ff 2 + ||w 2 || H 2 + H^IIhOIIwIIcc^tIjhi) ■ 
An integration over (0,r), Theorem 13.11 and Holder's inequality yield 

|(|| A^|| 2 + \\VifiW 2 ) + I J WMmdt < 27 7l ||w|| 2 2(0 T . H2) + fh||V||i a(0 , T; fl») 
+Ct|||(^,w,C)||| 2 + C , % /7||wi|| L 2 (0ir . ff 2 ) ||^|| 2 7([0iT];H2) 

+Cv / 7 : (||wi||i2( 0jT . J j2) + ||w 2 || L 2 (0ir;ff2) + \\iP 2 \\l 2 (0,t;H 3 )) \\w\\c(10,t];H1)- 

Hence, with a suitable choice of 771 , 772 and r, we have 

|(||A^|| 2 + HV^II 2 ) + I f H^ll^dt < L r |||(^,w,C)||| 2 , (57) 
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where < L T < 1. A comparison with (|55|) and Young's inequality yield 

IIA^ < C(\\<p t \\ 2 Hl + ||$!- $ 2 f Hl ), 
which, in view of (l26l) and (|5"7|) , guarantees 



imiw* < c / (ii^n^ + iia^ii^)^ < LriiKv.w.or. 

Jo 

From (05]) and (J3S) it follows that 

v t + vAv = Tr - T 2 , (58) 

where 

Ti — Ta = 7>[-kV • (V-01 (g) VV> + VV> <g> Vi/^) + AV'iV'itw + AV'iV , tw 2 + AtAV^w^ 
+Ai/'i(wi • Vi/'i)w + AV'(wi • VV'i)w2 + A'02(w • VV'i)w2 

+Al/>2 (W2 • V?/>) w 2 — (Vw)wi — (Vw"2)w]. 

Let us multiply (155)) by (v + Av) thus obtaining 

^[||v|| 2 + HlVvf] + ^(||Av|| 2 + ||Vv|| 2 ) < CUT! - T 2 || 2 . 

By means of the Sobolev embedding theorems and (|43l) . (PHI) inequalities, the same argu- 
ments used to prove (|53[) lead to the inequality 



||T X - T 2 || 2 < C[(||Vi||fl» + + VUWh + Utf + U^tMlVll^ 

+ l|w 2 || H 2 +i)||w||^i] + m\H\m +%Wfl» +%IIV'tllffi, 

for any 771 , 772, 773 > 0. An integration over (0, r) leads to 

i[||v|| 2 + H|Vv|| 2 ] + ^j\\\M\ 2 + ||Vv|| 2 )dt 

< CVr [II^iIU^O.t;^) + 11^2^2(0^.^3) 
+ ||'02t|U 2 (O,r ; H 1 )] \\4>\\c(0,t-M2) 

+Cy/f [||?Alt||L2(0,r;Hi) + 1^21^2(0^.^2)] HwH^o^.^i) 
+ GT [HV'IIcCO.t;^) + ll^llc(0,r;L2) + IM^o^i) 
+^l|| w llL2(0.r;H2) +?72||V'lli2(o^ ; H3) +773||^lli2 ( o !r;// i). 

Hence 

i[||v|| 2 + H|Vv|| 2 ] + ~ \\\Avf + HVv|| 2 )di < L T |||(^w,C)|||. (59) 



2 L J 2_ 
From ([37]) and gO), we obtain 

sfi t - 7 A/i = Ai - A 2) (60) 
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with 

Ai — A2 = —ipt — w i • V^> — w • V^2- 
We multiply (|60p by (/z — A/i) and we integrate over fi: 

|||| M ||l fl + 7 ||AH| 2 +7l|VMl| 2 

< ciihi 2 h1 +^(iiv t ii^ + \\^i\\umh + \\H\h\w\m)- 

Gronwall's inequality yields 



^Mm +7 y o IIA^Ir* < r?e° T [HVtlli^o^i) + l|wi||i 2(0 , r;H2 )|Mlc([o,T] ; H3) 
+IIV'2||| 2(0 , r; H3)l|w|| 2 7([0!r];H1) ] <i T |||(^w,C)||| 2 . (61) 

Therefore, collecting ([57)1 , (jSl?)) , (|6"Tj) and choosing r, 771,772,% small enough, we prove 
that 

IMIi^o.-nffS) + WftWh^gi) + \H\h(a,r-H 2 ) + IMIcao.T];^) < ^rllKV^w^)!!! 2 . 

The control on the remaining norms ||</9t||c([o,T];L 2 ) an d II v IIl 2 (o,t;L 2 ) is obtained by compar- 
ison with (|S"6")) and (|58p respectively. Therefore ([55)1 is proved. 

By means of a fixed point argument, the previous steps allow to prove that problem (P e ) 
admits a unique local solution {ip, v, /j) in X T , provided that r is small enough. □ 

Lemma 4.3 Any solution (<p £ ,v £ ,/i e ) of problem (P £ ) satisfies the estimates 

||v e W|| 2 + ||^(t)||^ +£ ||^W|| 2 

+ f [IIV^H 2 + ||v E ||^ + ||^|| 2 + W\\ 2 H2 ] dx < C Q (T), (62) 
Jo 

4f(t)\\h + I [\m 2 m+4vt\\ 2 ]dt<C n {T), (63) 



||A^(t)|| 2 + ||Vv £ (i)|| 2 

+ AllV^H 2 + 11^11^3 + IIAv^l 2 + K|| 2 ]di < C (T). (64) 
Jo 

for any t € (0, T), where Cq(T) is a positive constant depending on the inital data, the source 
f and the time T. 

Proof. By multiplying (|3"Tj) by ip\ + v e • V<p £ , we deduce the equality 

f ^l|V^|| 2 + (3\\cpl + v £ • V^|| 2 ~ «/V • V^)A^x 

+ / [(</) 3 + ^ £ + A(v £ )V-//](^+v £ - Vip e )dx = Q. (65) 
n 
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Now let us multiply (1521) by v e 

-A<^ e (<pf + v £ • V</)(v e ) 2 - f • v e ] dx = 0. 
Finally, we multiply (1331) by and integrate by parts, thus obtaining 



|Vv e || 2 + ( [kV ■ (V</ <g> \7ip e ) ■ v £ 
Jn 



e d 
2~dl' 



(66) 



(67) 



||^|| 2 + 7 ||V^|| 2 + f (ifil + v E • V^dx = 0. 

It is easy to show that since v e g _ff^ iv (Sl), the following identity holds 

/ V ■ (Vip e <Z)Vip E ) -v e dx = / (v £ • Vip £ )Aip £ dx. 
Jn Jn 

Therefore summing up equations (l6"5l) - (|67l) . we have 

/3||^+v e .V^|| 2 + i ,||Vv £ || 2 

+7l|V//|| 2 = - / [(^ £ ) 3 + u</]v e • \7tp £ dx+ I i-w £ dx. 
Jn Jn 



ld_ 
2~dl 



K ||V^|| 2 + ||v £ || 2 +s||^|| 2 + -||(^) 2 + W || 2 



(68) 



in Jn 
The first integral in the right hand side vanishes as a consequence of the identity 



f [& £ f + uLp £ ]v £ ■ \7tp £ dx = [ V 
Jn Jn 



E\4 



<f\2 



l dx 



and by applying the divergence theorem. Therefore, an integration of (|6"5|) over (0, t) and 
Poincare inequality provide 



«||V^(i)|| 2 + \W £ {t)\\ 2 +e\\lf{t)\\ 2 + \W{t? + u\\ 2 



/3||^+v-.V^|| 2 + -||Vv £ || 2 + 7 ||V^|| 2 



dt < C* , 



(69) 



where Co > depends on the norms of the initial data ||/Uo||> || v o||j Hyollif 1 and of the source 
||f||. In addition, an application of Young's inequality leads to 



ii^(i)ii 2 <c(ii^w 2 + w ir + i)<c , 



(70) 



for all te[0,T}. 

Now we multiply (l3"TT) by A(p £ . An integration by parts, Holder's and Young's inequalities 
imply 



e 1 1 2 r u n I ril,„e i Y7,„£||2 , I|v7,,e||2 , ||v7,„e||2 , ||/,„e\3 , „,,„e||2 



II A^ e || 2 ^ < c / \m + v £ ■ v^ii 2 + \\Vfi £ \\ 2 + \\v<p £ \\ 2 + \\(<p £ y + u<p £ \ 

o Jo 

+y £ (v £ n 2 ]dt. 



(71) 
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In view of Sobolev embedding theorem we obtain 

11(^)3 + ^e||2 + ||^ (v e )2|| 2 < CQ]^ + ||^||2 + ||^||^ || V E |||) (72) 

and by means of (|22|) . (l69l) and (j70|) , last term of (|72|) can be estimated as 

||^|| 2 ffI ||v £ ||| < C ||vl 8 / 3 ||v e ||^ 3 < C (||vl 8 + ||v^|| 2 ffl ) < C„(l + ||Vv E || 2 ). 
Therefore, on account of (|69|) - (l72l) . we have 

' \\f e \\ 2 H 2dt < C [ (H^ll 2 + \\A^\\ 2 )dt < C (T). 
Jo 

In addition the following inequality holds 

Mfdt < C [ (M+v E ■ V^ £ || 2 + ||v e ■ Vip e \\ 2 )dt 



< Co / (l + ||v e || 2 ||V^|| 2 )dt 

Jo 

< Co [ (l + \\v E \\ 2 m + \\V^\\ 2 H1 )dt<C . 

Jo 



Thus ([52]) is proved. 

By multiplying (f3"3"| in L 2 (il) by — A/+, we obtain 



||||V M E || 2 + 7||Am £ || 2 < jT + v £ • V^A/x^rr < CM + • V^|| 2 + J|A^|| 2 . 
Hence, 

e||V^(t)|| 2 +7 T W^ffdt < e||V Mo || 2 + C f M + v e • V^\\ 2 dt < C (T), (73) 



where last inequality follows from (I69|) . 
A comparison with (|33[) provides 

£ ||^|| 2 <C(||A^|| 2 + ||^+v £ .V^|| 2 ). (74) 

Finally we observe that (f3"Tj) implies 

||^H 2 < Co(||^|| 2 + 11^11^ + 11^11^). (75) 

Collecting (|73 > ([75 ]) . we get (53}. 

In order to prove ([54]) let us multiply (|3"Tj) by —Aip t , (15"2"j) by Av and integrate over f2. 
By means of Young's inequality we get 

~ N|A^|| 2 + llVv^l 2 ] + ^IV^II 2 + ^II^H 2 

< C [IIV^H 2 + ||(Vv £ )V^|| 2 + IRvv^KH 2 + ||(</) 2 V</|| 2 

+ ||V^|| 2 + ||V^(v e ) 2 || 2 + |b e (Vv e )v e || 2 + ||(VV</)V^|| 2 

+ ||^v e || 2 + ||^(v £ • V^K|| 2 + ||f|| 2 ] . (76) 
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Owing to (|2"2")l and (|62p , some of the terms of the right-hand side can be estimated as 

||(^) 2 V^|| 2 + ||V^|| 2 + ||V^) 2 || 2 < C (T)(||^||1, 2 + KUiO. (77) 
Similarly the remaining terms can be controlled as 



||^(v £ • V^KI^ < C||vl^T Hl ||V* < X(t)(h £ \\^ + \\<ff\\%>) (78) 
||^(Vv e )v e || 2 + ||(Vv E )V^|| 2 < CHVvlKll^ll 2 ^ l|v e || 2 + ||V^|| 2 ) 

<v\W e \\ 2 H2 + X (t)(\W £ \\ 2 m +l) (79) 
IKVV^KH 2 + ||(VV^)V^|| 2 < C||VV^|| 2 (||v e || 2 + ||V^|| 2 ) 

< ^11^11^3+ xW(l|v £ ||^i + 11^11^+1) (80) 

||^v-|| 2 < C\\^f H1 M\\l\W £ \\l < v\\ V^|| 2 + x(t)(h e \\m + 1), (81) 

where x is a L 1 function of time and r\ is a suitable positive constant. 

In order to evaluate the H 3 — norm of <p e , let us take the gradient of (j3"Tj) and obtain 

HVA^H < C [||V^|| + ||V//|| + ||Vv £ || 4 ||V^||4 + ||v £ || 4 ||VV^||4 
+lb e Hil|V^||e + l|V^|| + HV^IIeKH 2 ||Vv £ || 4 ||v £ || 6 ] . 

Interpolation inequalities (JUJ), (|2"21 and estimates ([52 ^) - ([55)l imply 

IIVA^H < m\\v E \\m + fj2||^||fl3 + C[||V^|| + ||V M £ ||] 
+C 



Hl ||^|| H2 + ||^||^ + l|v £ || ffl +||v e ||^ 3 



with rji,r/2 > 0. Choosing suitably r]2 and owing to (jl2o]) , (|62p . we prove 

Ibllla < C[||V^|| 2 + ||V^|| 2 ] + X (t)[l + \\v s \\ 2 H1 ] + m\W £ \\ 2 m, (82) 
Substitution of ([77 1) -([8"T ]) into (|75j!, use of {B2]), dZ3) and (JS2J provide 

~|N|A</|| 2 + HVv-ll 2 ] + ^IV^II 2 + P -\\Av*f + CM% 3 

< X (t)(||A^|| 2 + ||Vv £ || 2 + l). 

Gronwall's inequality and comparison with yield □ 

Proposition 4.1 Problem (P e ) admits at least a solution (ip £ , v e , 6 Xy. 

Proof. Theorem 14.11 guarantees existence of a solution (ip e ,v e , fi E ) defined in a small time 
interval (0,r). In order to extend this solution to the whole interval (0, T) we need the 
following uniform estimate of the solution 

\\<p e (i)\\ H 2 + \\v e (t)\\ m + \\n e (t)\\ m <K, te [0,T], (83) 

where if is a positive constant depending only on the global data ipo,Vo,fiQ and e, but 
independent of t. Inequality (l83|) is ensured by the estimates (|62)) - (l64l) of Lemma 14.31 
Therefore, by applying Theorem 14. 1[ after a finite number of steps we find a global solution 
of problem (P e ) in Xt- □ 
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5 Well posedness of the original system 
5.1 Existence of solutions 

Theorem 5.1 Let </? € H 2 (Vl), v € Hl iv (tt), fj, € H 1 ^). Then for any T > 0, there 
exists at least a solution (ip,~v,[i) £ Xt of prolem (P). 

Proof. Let (ip 5 , , v e ) € X T be a global solution of problem (P £ ). From the a priori 
estimates of Lemma 14.31 we deduce that 

(f £ is uniformly bounded in L 2 (0, T, H 3 (n)) n F x (0, T, H 1 ^)) 
v e is uniformly bounded in L 2 (0, T, # 2 (ft)) n ff x (0, T, £ 2 (0)) 
is uniformly bounded in L 2 (0, T, i? 2 (fi)) 
Ve^ 6 is uniformly bounded in H l (0, T, L 2 (fi)). 

As a consequence there exists a subsequence, denoted also (ip £ , v e , fi e ) such that 



as e — > 0. In particular, as a consequence of Aubin's theorem, we have 

(p £ -> tp strongly in L 2 (0, T, H 2 (fl)) 
v £ ^v strongly in ^((LT,^ 1 ^)). 

The previous convergences allow us to pass to the limit as e — > into (|3"Tj) - (l3"3"|) and to obtain 



5.2 Uniqueness of solution 

Theorem 5.2 Problem (P) admits a unique solution (ip,v,fi) G Xt- 

Proof. First we notice that by repeating the same arguments of Lemma 14.31 with e = 0, one 
can easily show that any solution (c/2, v, (i) of problem (P) satisfies the following estimate 

\W(t)\\m + \\<p(t)\\h + f [IMIk + Wftfm] dt < Co(T). (84) 
Jo 

Let (ipx, vi, /ii), (v2,V2,/i2) € At be two solutions to problem (P) with the same initial 
data and source f . We consider the differences 



f e -> f 
ft -> P< 



vf -> v t 
e/4 -> 



weakly in L 2 (0,T, # 3 (ft)) 
weakly in L 2 (0, T, F x (0)) 
weakly in L 2 (0,T, # 2 (ft)) 
weakly in L 2 (0,T, L 2 (fi)) 
weakly in L 2 (0,T, iJ 2 (f2)) 
strongly in L 2 (0, T, L 2 (f2)) 



(EHj-daoj). 



□ 



v = vi — v 2 , 



22 



satisfying the following problem 



with 



ipt = — v • V</?i — V2 • V(f + yA/U (85) 
v t = —uAv — V[kV ■ (V</?i ® V</? + V</5 V^ 2 ) + \<pi((fu + vi ■ V^i)vi 
-\(f2(^2t + v 2 ■ V(^ 2 )v2 + (Vv)vi + (Vva)v] (86) 



H = (3ip t — nAip + /3vi • + /3v • Vip 2 + (<Pi + V? 2 + Vi^f + wp + \ipv\ 
+Av? 2 (vi + v 2 ) • v. (87) 



We append to ([80 ]) - ([87 ]) the initial conditions 

y>(x,0)=0, v(x,0) = 0, a.e. x G O. 
We multiply equation (|85[) by (p and we integrate over O. Thus, we obtain 
1 d f 

2j t M 2 =-J t( v ' V ^ + V2 ' V ^ + ' dx - 
By means of the Holder and Young inequalities and in view of ([84]) , we have 

^IMI 2 <l|v|| 2 + |||v^| 2 + c|M|^. 

Now let us multiply (|55| in L 2 (f2) by [i, thus obtaining 

/ [(^ t /i + (vV( y 3 1 +v 2 -V^)M + 7|V ( u| 2 ]dx = 0. (89) 
The first term of the integral may be rewritten by substituting expression (1571) as 



/ 



tptVdx = ~\\Vtp\\ 2 + f3y t \\ 2 + I [/3(v • V^i + v 2 • Vtp) + ^{ifj + + tpitpa) 



n 2 dt jfi 

+u<p + Xfvf + A</? 2 (vi + v 2 ) ■ v]ip t dx. 

A substitution into ([89]). Holder's inequality and a priori estimate (|84|) yield 

||||V^|| 2 +/3||^|| 2 + 7l|V/,|| 2 

< C[(l + HvallooJII^IlHi + ||v|| 4 ] IbtH + C[||v|| 4 + I|v 2 |UI|V^||]||m||. 
Young's inequality, (f2T|) and ([23]) lead to the estimate 

f|l|V^|| 2 + ^||^|| 2 + g7l|VH| 2 (90) 
< C(l + Hvalll,)^!!^ + q|v|| 2 + T^||Vv|| 2 + m M\ 



where f]\ is a suitable positive constant. 
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By multiplying (|87[1 by /i and accounting for ([84")h we infer that 

ll^ll 2 < CUVmIIIIV^II + C[\\<pt\\ + || Vl |MMI^ + l|v||ifi + IMMIMI- 

H6I1C6 

IIHI 2 < c[||V/i|| 2 + ll^f + IM^lMlk + Hi 2 ,, + Mir}. (91) 

By multiplying (|86l) by v we obtain 

1 d 



where 



J|v|| 2 + H|Vv|| 2 < h + h + h + h, (92) 

h = k I (V^i ® V<y9 + Vip ® Vip 2 ) ■ Vvefe 
Jn 

,/n 

I 3 = X [<p(vi ■ Vvi)v a + y> 2 (vi • Vyi)v + <^ 2 (v ■ V^i + v 2 ■ V^)v 2 ] • vdx 
Jn 

h = ~ [(Vv)vi + (Vv 2 )v] • vdx. 
Jn 

Holder's, Young's inequalities, (|2"U|). (|2"T1) and a priori estimate allow us to estimate Ji, 
namely 

/i < «(I|V^|| 4 + ||V^|U)||V^|| 4 ||Vv|| < ^||Vv|| 2 + V2 y\\h + C||V^|| 2 . (93) 
Similarly, we obtain 

h < A[||<y5||//i||(piJ||vi|| H i + ||^ 2 || H i||^ t ||||vi|| H i + ||¥?2||ffi||v?2t||||v|| 4 ]||v||4 



<^l|Vvf + |||^|| 2 + C(||vf + 11^) (94) 



and 



h < \[\\<p\\# living \\ipi W& + (||vi||hi + ||v2|U0ll^||Hi||^i|kH|v||4 
+||^2||hi||v2|| 2 ji||V93|| 4 ]||v||4 

< ^l|Vv|| 2 + m M 2 H2 + C(M\m + llvll 2 ). (95) 
Finally, last integral can be controlled as 

h < llVvllllvilUHvlU + ||Vv 2 ||||v|| 2 < ^||Vv|| 2 + C(|| Vl || 2 ffl + ||v 2 || 2 f2 )||v|| 2 . (96) 
From (JSTJ) it follows that 

||A<p|| < C \\fi\\ + |M|jji|l¥>i + ipnp 2 + + + |M|jHI v iIIjp 

+ ||¥>2||ifl||vi + v 2 || H i||v|| 4 + \\ip t \\ + ||vi|| ff 2||V^|| + ||V^ 2 || 4 ||v|| 4 

< C(\\fx\\ + \\ip\\ H i + \\<pt\\ + ||v|| 4 ) + C||v 2 || H 2||V^||. 
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As a consequence, by means of (l9"T1) . we have 

IMI^<c(|M| 2 + ||A^|| 2 ) 

< C [||V/*|| 2 + (|| Vl ||| 2 + ||v 2 ||| 2 + 1)11^111, +||^|| 2 + \\A 2 m] ■ (97) 

Adding inequalities and (j9"2")l . accounting for ([^T]). (|^7|) and choosing 771,772 small 

enough, we prove the estimate 

\j t [\M 2 + kIIV^I 2 + ||v|| 2 ] + ^i^il 2 + Jiv^ii 2 + pw|| 2 < h(t)(\\<P\\h + HI 2 ) 

where 

^)=C(l|vi|| 2 ff2 + ||v 2 || 2 ff2 + l) 
is a L 1 — function of time. Thus, Gronwall's inequality proves 

ip = 0, v = 0. 

Accordingly, from (|87|) it follows that /x = and we reach the conclusion. □ 
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